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Improved  Sand  Transport  Model 
for  LTFATE  (Version  2.0) 


PURPOSE:  Long-Tenn  FATE  of  dredged  material  (LTFATE)  model  is  a  combined  local 
hydrodynamic  and  sediment  transport  model  used  to  determine  the  long-  and  short -lerm  stability 
of  dredged  material  mounds.  This  Technical  Note  (TN)  describes  a  more  accurate,  replacement, 
sand  transport  submodel  that  has  been  incorporated  into  LTFATE. 

BACKGROUND:  The  LTFATE  model  includes  two  sediment  transport  methods,  cohesive  and 
noncohesive.  The  user  specifies  which  type  of  transport  is  being  simulated  in  the  input  file.  This 
decision  is  based  on  an  understanding  of  the  dredged  material  mound  composition  being  analyzed. 
Predominately  sandy  mounds  where  interparticle  forces  between  grains  are  small  comparec.  to 
gravitational  forces  are  modeled  as  noncohesive.  Sediments  that  contain  sufficient  fines  are  modeled 
as  cohesive.  A  description  of  the  LTFATE  cohesive  sediment  transport  routines  is  provided  in 
Gailani  (1998). 

The  amouni  of  fine-grain  material  needed  to  make  sediment  behave  in  a  cohesive  manner  is  a 
function  of  several  variables,  including  mineralogy,  organic  content,  and  grain-size  distribuEion. 
There  is  no  formula  available  for  making  this  decision.  However,  for  screening  level  efforts,  it  can 
be  assumed  that  sediment  containing  less  than  approximately  five  to  eight  percent  (by  weight) 
smaller  th  in  10  pm  and  a  majority  of  the  material  greater  than  80  pm  (sand-size)  can  be  considered 
noncohesive.  This  does  net  hold  true  in  all  cases.  For  some  very  cohesive  fine  sediments,  such  as 
bentonite,  fractions  less  than  two  percent  can  make  a  predominately  sand  ( 1 00  pm)  .sediment  behave 
as  a  cohesive  material  (Gailani  et  al.  200 1 ).  Erosion  experiments  are  the  best  method  for  determining 
if  the  sediment  of  interest  is  cohesive.  However,  general  knowledge  of  the  sediment,  visual  analysis, 
or  rules  of  thumb  can  be  used  when  erosion  data  are  not  available.  LTFATE  Version  2.0  requires 
that  the  user  has  determined  a  priori  if  the  sediments  are  cohesive  or  noncohesive . 

LTFATE  Version  1 .0  (Scheffner  et  al.  1 995;  Scheffner  1 996)  simulated  sand  transport  using  a  total 
load  formula  developed  by  Ackers  and  White  (1973).  The  model  simulates  this  total  load  movement 
as  bed  load.  The  version  of  LTFATE  described  in  this  TN  includes  methods  for  estimating  both 
suspended  and  bed  load  of  sandy  sediments.  The  method  developed  by  van  Rijn  (van  Rijnl989a,b; 
1993)  has  been  incorporated  into  the  LTFATE  model.  Compared  to  the  Ackers  and  Wliite  method, 
the  van  Rijn  method  has  been  more  rigorously  tested  and  compared  to  data  in  high-energy  and 
wave-dominated  regimes.  Therefore,  Version  2.0  permits  more  accurate  simulation  of  sand  trans¬ 
port,  particularly  during  high-energy  events,  where  some  sand  is  in  suspension  and  can  be  carried 
long  distances  before  redeposition. 

Van  Rijn  Sand  Transport  Algorithms.  The  van  Rijn  sand  transport  routines  incorporated  into 
LTFATE  Version  2.0  were  designed  and  calibrated  for  nonbreaking  wave  conditions.  Therefore, 
they  do  not  include  suspension  and  transport  induced  by  breaking  waves  and  are  not  appropriate 
for  the  surf  zone.  The  routines  do  estimate  wave  asymmetry  and  are  therefore  designed  for  nearshore, 
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nonbre^ng  wavj&  cpnditipns.  The  method  developed  by  van  Rijn  (1993;  1989a,b)  is  time-averaged 
over  a  wave  cycle  for  Suspended  load  and  time  varying  for  bed  load.  Wikramanayake  and  Madsen 
(1994)  and  otlicm  criticize  tim  averting  as  missing  the  spikes  in  concentration  caused  by  passage 
of  the  trough  jand  pe^k  of  die  LTFATE  grids  have  minimum  cell  size  on  the  order 

of  tens  of  ipeters  and  -time-steps  of  a  wave  period  or  larger.  Breakdown  into  smaller  grid  cells  or 
time-steps  that  would-be  required  to  monitor;  concentration  variation  over  a  wave  cycle  is  not 
practical.  Im  additioti;  the  ^opted  method  has  been  shown  to  reasonably  simulate  suspended 
transport  and  timeiaveraged  suSp^ded  concentrations  (van  Rijn  1993). 

The  method,  as  coded  ih  the  van  Wjn  TRANSPOR  routine  for  transport  at  a  point  (van  Rijn  1 993), 
requires  as  inpiitV  water  depth,  depth-averaged  ciirrent  velocity,  grain-size  distribution,  angle 
between  wave  and  current  direction,  cturent  and  wave  related  bed  roughness,  fluid  temperatm  e,  and 
salinity.  The  niodM  willicalculate  wave  asymmetry  based  on  wave  and  water  depth  conditions.  Wave 
asymmetry  is  assuih^  to  afiect  bed  load,;  but  not  suspended  load  horizontal  fluxes. 

Dimensionless  Pai:ameters.  Three  dimensionless  parameters  are  used  in  the  van  Rijn  method. 
The  dimensionleSsipaitiC,l5  parameter  Zj*  reflects  the  influence  of  gravity,  density,  and  viscosity  on 
sediment  particles^ 

D*=45o|^(s-l^g/v^]  ^  ■  (1) 

Where  </5o  is  the  sediment  bed  median  particle  diameter,  s  is  the  specific  density  of  the  sediment,  g 
is  the  acceleration  pf  gi^Yity,  and  v  is  the  kinematic  viscosity  of  water.  The  Shields  parameter  for 
initiation  of  movement  is  used  to  calculate  wave  and  current  related  critical  shear  stresses  and  can 
be  estimated  from  grain  size: 

t-'  -■  ■  Y  ''  '  ■  :  '  \ 

Qcr-bDS  I"  ’  :  '  (2) 

where  b  and  c  are  Coefficients  based  on  the  value  ofD*  such  that  the  value  of  the  Shields  parameter 
approximate  toe  weU-kno^^^^  curve  (van  Rijn  1993).  The  dimensionless  bed  shear  stress 

parameter  for  bed^oad  ti^port  T,  and  at  height  a  above  toe  bottom  (a  is  toe  height  of  toe  bed 
load  layer)  are  estimated  respectively  by: 


where  is  toe  wave-current  interaction  coefficient  which  is  a  function  of  bed  roughness,  wave 
boimdary  layer  thickness,  and  water  depth;  jXj.  is  toe  current  efficiency  factor,  which  is  a  function 
of  bed  rouglmess  arid  grain-sizevdisttibution;  and  ^  are  the  wave  efficiency  factors  at  toe  bed 
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and  at  height  a  above  the  bed,  respectively;  |\,  and  ^  ^  function  of  the  grain-size  distribution, 

bed  roughness  and  wave  conditions.  The  methods  for  estimating  the  four  parameters  are  complex 
and  the  user  is  referred  to  van  Rijn  (1993)  for  additional  details.  In  addition,  is  the  current  related 
bed  shear  stress;  ^  is  the  wave  related  bed  shear  stress;  and  x^  j,,.  is  the  bed  critical  shear  stress  for 
the  initiation  of  sediment  resuspension.  These  shear  stresses  are  calculated  in  the  code  and  are  based 
on  currents,  wave  conditions,  and  friction  factors.  Other  parameters  calculated  in  the  code  include; 
backward  and  forward  orbital  velocities,  active  and  boundary  layer  thickness,  friction  factors,  and 
sediment  mixing  coefficients.  Again,  the  user  is  referred  to  van  Rijn  (1989a,b;  1993)  for  details  on 
their  derivation. 

Sediment  Fluxes.  The  described  parameters  are  then  used  to  estimate  time-varying  bed  load  and 
time-averaged  suspended  load  fluxes.  These  fluxes  are  then  input  to  the  LTFATE  model  to  move 
sediment  from  cell  to  cell.  The  time-varying  bed-load  transport  rate  (m^/m/s)  is  estimated  as  (van 
Rijn  1993): 
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where  7=1-  h  is  the  water  depth  (m),  is  the  significant  wave  height  (m),  p  is  water 

density,  is  the  sediment  density,  is  the  critical  bed-shear  stress  according  to  Shields  (N/m^), 
and  x'ly  cw  is  the  grain-related  instantaneous  (time-varying)  bed  shear  stress  due  to  currents  and 
waves  (N/m^).  In  addition,'ffc,cH'  is  related  to  both  x^  ^  and  x^  but  is  not  the  sum  of  the  two  shear 
stresses.  Summing  wave  and  current  related  shear  stresses  has  been  shown  to  underestimate  total 
bottom  shear  stress  (Grant  and  Madsen  1979).  The  method  used  by  van  Rijn  accounts  for  current 
and  wave  processes  interactions  to  create  a  higher  shear  stress  than  the  two  components  separately. 
Time-averaged  values  of  qi,(t)  are  obtained  by  averaging  over  a  wave  period. 

Once  it  has  been  determined  that  the  bottom  shear  stress  exceeds  the  critical  values  for  suspension, 
then  an  equilibrium  volumetric  reference  concentration  at  a  reference  height  z = a  above  the  sediment 
bed  is  calculated  as: 


y^a,eq  ~  9-915 
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In  addition,  the  concentration  profile  in  the  water  column  is  calculated  by  ': 


(6) 


_ [l-c(z)]  ^c(z)  vv, _ 

e,.,,.„,(z)  |l  +  [c(z)/co]  -2[c(z)/c„]  "'"‘j 

^  (7) 

where  c  is  the  volume  concentration,  Cg  is  the  maximum  volume  concentration  (=  0.65),  is  the 
particle  settling  speed,  and  is  the  combined  current/wave  sediment  mixing  coefficient  (van 
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Rijn  1993).  8c, jw  is  a;<?omplicate4  function  of  wave,  current,  water,  and  sediment  conditions.  The 
user  is  referred  to  van  Rijn  (1993)  for  further  details.  The  concentration  profile  is  critical  in 
determining  die  aihoxint  bf  suspended  material  that  settles  in  each  time-step.  Deposition  is  based  on 
the  near-bottom  OoncentratiOp  profile,  not  the  vertically  averaged  concentration.  Settling  speed, 
(cm/s),  is  estimated  as  a  function  of  particle  size  by  (Cheng  1997): 
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Bed  load  apd  suspended  tiuxes  across  cell  boundaries  are  then  calculated  and  concentration  and 
sediment  bed  elev^on'chaiige  updated.  Suspended  load  flux  in  the  direction  of  flow  is  calculated 
by:.  ^ 

= 

^  ...a  ;/  V  ■ 

where  « is  the  vertic^ly  averaged  velocity  determmed  in  the  hydrodynamic  section  of  LTFATE. 
Suspended  and,  bpd-lOad  fiukeS^i^^  divided  into  x  and  y  direction  fluxes  within  LTFATE.  After  bed 
load  and  stispended  fluxes  across  cell  boundaries  are  calculated,  concentration  and  sediment  bed 
elevation  ohange  are  updated.  ThiS  provides  a  new  value  of  Q.  The  resuspension  flux  fi'om  the  bed 
to  the  water  cdlutnn  at  the  next  time-step  is  calculated  using  (van  Rijn  1993): 


(10) 

Where is  the' net  resuspension  flinc  Ifrom  the  sediment  bed  and  is  the  suspended  sediment 
concentration  a:t  z  -  4,  calculated  by  the  fluxes  across  the  cell  elements.  If  Q  <  Q  additional 
sediment  can  be  carried/by  ffew,  so  that  erosion  occurs  and  Ey^  >  0.  Conversely,  if  the 

carrying  capacity  of  the  water  colliinn  at  a  particular  bed  shear  velocity  has  been  exceeded,  which 
means  that  <  P,anti  deposition  Can  occur,  even  though  the  critical  bed  shear  velocity  has  been 
exceeded.  Wheii  l^d  shear  velocity  is  less  than  the  critical  value  or  the  resuspension  rate 

is  zero  and  noncohesive  sediments  in  the  water  column  can  be  deposited  on  the  sediment  bed.  The 
deposition  flux  is  calculated  as  Q  multiplied  by  Wj. 


RESULTS  01®  LTFATE  MODEL:  A  sample  simulation  was  performed  and  results  shown  here 
to  demonstrate  the  new  LTFATE  sand  code.  Cell  size  was  set  at  15.24  m  (50  ft)  for  this  simulation 
and  the  grid  was  comprised  of  61x61  cells  in  the  x  and  y  directions.  Water  depth  surrounding  the 
mound  was  15.24;#  (50  ft).  Storm  wave  and  current  conditions  are  shown  in  Figure  1 .  U  velocity 
is  in  the  x-direclioh.  Becausf  this  was  a  storm  simulation,  it  was  considered  “short  term”  and  data 
were  input  to  the  model  throng  a  “storm.dat”  file.  The  user  is  referred  to  Scheffiier  et  al.  (1995) 
for  formatting  of  this  in^t  file.  A  profile  of  the  mound  cross  section  in  the  x-direction  before  and 
after  the  storm  event  is  shown  in  Figure  2  The  combination  of  large  waves  and  high  current  in  the 
negative  x-direction  during  this  storm  resulted  in  migration  of  the  mormd  in  the  negative  x-direction. 
In  addition,  mounti  peaks  near  the  leading  edge  were  reduced  in  height. 
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